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2i ■ Abstract 

-)— » 

^j I I review explanations for the three neutrino anomahes (solar, atmospheric 

[■^^ I and LSND) which go beyond the "conventional" neutrino oscillations induced 

by mass-mixing. Several of these require non-zero neutrino masses as well. 

> 

t:j- . I. INTRODUCTION 

(N 

As is well-known, it is not possible to account for all three neutrino anomalies with just 

0> ■ the three known neutrinos (z/g, z/^ and Ut-). If one or more of them can be explained in some 

other way, then no extra sterile neutrinos need be invoked. This is one motivation for exotic 

scenarios. In any case it is important to rule out all explanations other than oscillations in 

order to establish neutrino mixing and oscillations as the unique explanation for the three 



D I observed neutrino anomalies. 

"^ ' I should mention that, in general, some (but not all) non-oscillatory explanations of the 

neutrino anomalies will involve non-zero neutrino masses and mixings. Therefore, they tend 
S^ . to be neither elegant nor economical. But the main issue here is whether we can establish 

oscillations unequivocally and uniquely as the cause for the observed anomalies. 

I first summarize some of the exotic scenarios and then consider each anomaly in turn. 

II. MIXING AND OSCILLATIONS OF MASSLESS u'S 

There are three different ways that massless neutrinos may mix and even oscillate. These 
are as follows. 

1. If flavor states are mixtures of massless as well massive states, then when they are 
produced in reactions with Q-values smaller than the massive state; the flavor states are 
massless but not orthogonal |I[|. For example, if u^ = Yli=iUeii^i, i^fj, = Z]i=i ^/^i^i and 
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mj = for i = 1 to 3 but m4 = 50 GeV; then 'Vg" produced in /3-decay and 'V^" produced 
in TT-decay are niassless but not orthogonal and 

(^e I ^m) = -U:, U^ (1) 

On the other hand Ue and u^^ produced in W decay will be not niassless and will be 
more nearly orthogonal. Hence, the definition of fiavor eigenstate is energy and reaction 
dependent and not fundamental. Current limits on orthogonality of u^, v^ and v^. make 
it impossible for this to play any role in the current neutrino anomalies. 

2. When Flavor Changing Neutral Currents as well as Non- Universal Neutral Current 
Couplings of neutrinos exist, for example: 

AG 

+ 4"?^ {^^L ^/^ ^^L - ^rLlt, i^tl} {ql 7/. Ql} (2) 

Such couplings arise in R-parity violating supersymmetric theories 0]. In general, in these 
theories, one expects neutrino masses at some level here it is assumed that the FCNC 
provides the dominant effect. In this case the propagation of I'e and Ur in matter is described 
by the equation: 

(3) 

where a = GpNe.lS = egGpNg and 7 = e'^GpNg. There is a resonance at a + 2/5 = or 
eg = —^Nf./Ng and u^ can convert to u^ completely. This matter effect which effectively mixes 
flavors in absence of masses was flrst pointed out by Wolfenstein in the same paper IQ where 
matter effects were flrst discussed. The resonant conversion is just like in the conventional 
MSW effect, except that there is no energy dependence and u and u are affected the same 
way. This possibility has been discussed in connection with both solar and atmospheric 
anomalies. 

3. (a) Flavor violating Gravity wherein it is proposed that gravitational couplings of 
neutrinos are flavor non-diagonal [^ and equivalence principle is violated. For example, ui 
and 1^2 may couple to gravity with different strengths: 

Hgr = fiGE<j) + f2GE(j) (4) 

where (p is the gravitational potential. Then if z/g and z/^ are mixtures of z/i and 1^2 with a 
mixing angle 6, oscillations will occur with a flavor survival probability 

P = 1 - sin^2dsin'^{-5f(pEL) (5) 

when (f) is constant over the distance L and Sf = fi — /2 is the small deviation from 
universality of gravitational coupling. Equivalence principle is also violated. 

(b) Another possibility is violation of Lorentz invariance wherein all particles have 
their own maximum attainable velocities (MAV) which are all different and in general also 





different from speed of light: then if z/i and V2 are MAV eigenstates with MAV's vi and V2 
and Vf, and v^ are mixtures of vi and z/2 with mixing angle 6*, the survival probability of a 
given flavor is 

P = 1 - sin^2e sin^ (-5vEL\ (6) 

where 5v = Vi — V2- 

As far as neutrino oscillations are concerned, these two cases are identical in their de- 
pendence on LE instead of L/E as in the conventional oscillations |^. 



III. NEUTRINO DECAY 

Neutrino decay |^ implies a non-zero mass difference between two neutrino states and 
thus, in general, mixing as well. I will consider here only non-radiative decays. We assume 
a component of fa, i-e., ^'2, to be the only unstable state, with a rest-frame lifetime tq, and 
we assume two flavor mixing, for simplicity: 

z/^ = cos9h'2 + sinOvi (7) 

with m2 > mi. From Eq. (2) with an unstable z/2, the z/q, survival probability is 

P„„ = sin^e + cos^^exp(-aL/P) (8) 

+ 2sin^ecos^ee^^{-aL/2E)cos{6m^L/2E), 

where 5m? = m\ — m\ and a = 1112 /tq. Since we are attempting to explain neutrino data 
without oscillations there are two appropriate limits of interest. One is when the 5m? is so 
large that the cosine term averages to 0. Then the survival probability becomes 

p^^ = sin^e + cos^^exp(-aL/P) (9) 

Let this be called decay scenario A. The other possibility is when 5ir? is so small that the 
cosine term is 1, leading to a survival probability of 

p^^ = (sin^e + cos'^eexp{-aL/2E)f (10) 

corresponding to decay scenario B. Decay models for both kinds of scenarios can be con- 
structed; although they require fine tuning and are not particularly elegant. 

IV. LSND 

In the LSND experiment, what is observed is the following Pj. In the decay at rest 
(DAR) which is /i"*" — > e'^UgUf^ which should give a pure z/g signal, there is a flux of z/^s at 
a level of about 3.10"^ of the u'^s. (There is a similar signal for u'^s accompanying u'^s in 
the decay in flight of /i~ — > e~z/eZ/^). Now this could be accounted for without oscillations 
||T0| provided that the conventional decay mode /x"*" -^ e"'"z/eZ/^(yU~ -^ e^Uf,v^) is accompanied 
by the rare mode /i^ -^ e"^z/eX (yU~ -^ e^z/gX) at a level of branching fraction of 3.10"^. 



Assuming X to be a single particle, what can X be? It is straight forward to rule out X as 
being (i) u^ (too large a rate for Muonium-Antimuonium transition rate), (ii) z/g (too large 
a rate for FCNC decays of Z such as Z ^ fie + fie) and (iii) u^- (too large a rate for FCNC 
decays of r such as r — >• /lee). 

The remaining possibilities for X are z/„ or z/gteriie- No simple models exist which lead to 
such decays. Rather Baroque models can be constructed which involve a large number of 
new particles |TT]]. 

Experimental tests to distinguish this rare decay possibility from the conventional oscil- 
lation explanation are easy to state. In the rare decay case, the rate is constant and shows 
no dependence on L or E. 

One can also ask whether the Ue events seen in LSND could have been caused by new 
physics at the detector; for example a small rate for the "forbidden" reaction z/^+p — > n+e"*". 
This is ruled out, since a fractional rate of 3.10"^ for this reaction, leads (via crossing) to a 
rate for the decay mode ii'^ -^ e'^i^n in excess of known bounds. 

V. SOLAR NEUTRINOS 

Oscillations of massless neutrinos via Flavor Changing Neutral Currents (FCNC) and 
Non Universal Neutral Currents (NUNC) in matter have been considered P,p!2| as explana- 
tion for the solar neutrino observations, most recently by Babu, Grossman and Krastev |]13| . 



Using the most recent data from Homestake, SAGE, GALLEX and Super-Kamiokande, they 
find good fits with e„ ~ lO'^ and e'^ ~ 0.43 or e^ ~ 0.1 to 0.01 with e'^ ~ 0.57. 

Since the matter effect in this case is energy independent, the explanation for different 
suppression for ^B, ^ Be and pp neutrinos as inferred is interesting. The differing suppression 
arises from the fact the production region for each of these neutrinos differ in electron and 
nuclear densities. This solution resembles the large angle MSW solution with the difference 
that the day-night effect is energy- independent. 

The massless neutrino oscillations described in Eq.(6) also offer a possible solution to 
solar neutrino observations. There are three solutions. Two are characterized by sin^ 26 ~ 
2.10"^, 5v/2 ~ 6.10"^^ and sin^ 20 ~ 0.7, 5v/2 ~ 10"^^; these are analogs of the small angle 
MSW and large angle MSW solutions fl^. These are now ruled out [|l5l (for both Ve — ^fj, 
and z/e - Ur) by the recent NUTEV data p[. The third one at sin^ 20 ~ 1 and |^ ~ lO'^^ 



(the analog of the vacuum solution) is the only one allowed |17]. This possibility can be 
tested in Long Baseline experiments. 

The possibility of solar neutrinos decaying to explain the discrepancy is a very old sug- 
gestion []TB| . The most recent analysis of the current solar neutrino data finds that no good 



fit can be found: Uei ~ 0.6 and Ti, (E=10 MeV) ~ 6 to 27 sec. come closest [|T^. The fits 



become acceptable only if the suppression of the solar neutrinos is energy independent as 
proposed by several authors [^ (which is possible if the Homestake data are excluded from 
the fit). The above conclusions are valid for both the decay scenarios A as well as B. 



VI. ATMOSPHERIC NEUTRINOS 

The massless FCNC scenario has been recently considered for the atmospheric neutrinos 
by Gonzales-Garcia et al. |^ with the matter effect supphed by the earth. Good fits 
were found for the partially contained and multi-GeV events with e^ ~ 1, e'^ ~ 0.02 as 
well eg ~ 0.08, e' ~ 0.07. The fit is poorer when higher energy events corresponding to 



up-coming muons are included p2|j23[| . The expectations for future LBL experiments are 
quite distinctive: for MINOS, one expects P^^ ~ 0.1 and P^^ ~ 0.9. The story for massless 
neutrinos oscillating via violation of Equivalence Principle or Lorentz Invariance is very 
similar. Assuming large z/^ — u^ mixing, and 6v/2 ~ 2.10~^^, a good fit to the contained 



events can be obtained [|2J]; but as soon as multi-GeV and thrugoing muon events are 



included, the fit is quite poor P^ 

Turning to neutrino decay scenario A, it was found that it is possible to choose 6 and a to 
provide a good fit to the Super-Kamiokande L/E distributions of u^^ events and z/^/z/g event 
ratio [^]. The best-fit values of the two parameters are cos'^9 ~ 0.87 and a ~ IGcV/De, 
where De = 12800 km is the diameter of the Earth. This best-fit a value corresponds to a 
rest-frame 1^2 lifetime of 

To = m2/a ~ j^ X 10-i°s. (11) 

However, it was then shown that the fit to the higher energy events in Super-K (especially 



the upcoming muons) is quite poor p2| , p3 



In all these three cases, the reason that the inclusion of high energy upcoming muon 
events makes the fits poorer is very simple. The upcoming muons come from much higher 
energy u' s and although there is some suppression, it is less than what is observed for lower 
energy events at the same L (zenith angle). This is in accordance with expectations from 
conventional oscillations. The energy dependence in the above three scenarios is different 
and fails to account for the data. In the FCNC case there is no energy dependence and 
so the high energy u' s should have been equally depleted, in the FV Gravity (or Lorentz 
invariance violation) at high energies the oscillations should average out to give uniform 50% 
suppression and in the decay A scenario due to time dilation the decay is suppressed and 
there is hardly any depletion of u'^s. 

Turning to decay scenario B, consider the following possibility [Q. The three weak cou- 
pling states v^iVti ^s (where Ug is a sterile neutrino) may be related to the mass eigenstates 
Z/2Z/3Z/4 by the approximate mixing matrix. 

/ cos 9 sin 6* \ / z/2 \ 




— sin 6 cos 6 
V 1 



^3 



(12) 
and the decay is z/2 — > z/4 + J. The electron neutrino, which we identify with z/i, cannot mix 



very much with the other three because of the more stringent bounds on its couplings |^ 
and thus our preferred solution for solar neutrinos would be small angle matter oscillations. 
Then the 6m23 in Eq. (1) is not related to the 6m24 in the decay, and can be very small, 
say < 10^^ eV^ (to ensure that oscillations play no role in the atmospheric neutrinos). In 
that case, the oscillating term is 1 and -P(z/^ —>■ v^) becomes 



P[u^ -^ u^) = (sin^ e + cos^ ^e-"^/2£;)2 (13) 

This is identical to Eq. (13) in Ref. |^. 

In order to compare the predictions of this model with the standard v^ *-^ Vr oscillation 
model, we have calculated with Monte Carlo methods the event rates for contained, semi- 
contained and upward-going (passing and stopping) muons in the Super-K detector, in the 
absence of 'new physics', and modifying the muon neutrino flux according to the decay or 
oscillation probabilities discussed above. We have then compared our predictions with the 
SuperK data p5[, calculating a x^ to quantify the agreement (or disagreement) between 



data and calculations. In performing our fit (see Ref. [^ for details) we do not take into 



account any systematic uncertainty, but we allow the absolute flux normalization to vary as 
a free parameter (3. 

The 'no new physics model' gives a very poor fit to the data with x^ = 281 for 34 d.o.f. 
(35 bins and one free parameter, f3). For the standard z/^ <-^ z/^ oscillation scenario the best fit 
has x^ = 33.3 (32 d.o.f.) and the values of the relevant parameters are Am^ = 3.2 x 10~^ eV^, 
sin^ 29 = 1 and /? = 1.15. This result is in good agreement with the detailed fit performed by 
the SuperK collaboration [^ giving us confidence that our simplified treatment of detector 
acceptances and systematic uncertainties is reasonable. The decay model of Equations (3) 
and (4) above gives an equally good fit with a minimum x^ = 33.7 (32 d.o.f.) for the choice 
of parameters 

T^/m^ = 63 km/GeV, cos^ 9 = 0.30 (14) 

and normalization jS = 1.17. 

In Fig. 1 we compare the best fits of the two models considered (oscillations and decay) 
with the SuperK data. In the figure we show (as data points with statistical error bars) the 
ratios between the SuperK data and the Monte Carlo predictions calculated in the absence 
of oscillations or other form of 'new physics' beyond the standard model. In the six panels we 
show separately the data on e-like and /i-like events in the sub-GeV and multi-GeV samples, 
and on stopping and passing upward-going muon events. The solid (dashed) histograms 
correspond to the best fits for the decay model (z/^ ^^ v^ oscillations). One can see that the 
best fits of the two models are of comparable quality. The reason for the similarity of the 
results obtained in the two models can be understood by looking at Fig. 2, where we show 
the survival probability -P(t'^ -^ v^) of muon neutrinos as a function of LjEy for the two 
models using the best fit parameters. In the case of the neutrino decay model (thick curve) 
the probability -P(z/^ — »■ Vf^ monotonically decreases from unity to an asymptotic value 
sin^ 9 ~ 0.49. In the case of oscillations the probability has a sinusoidal behaviour in LjE^. 
The two functional forms seem very different; however, taking into account the resolution 
in LjEy^ the two forms are hardly distinguishable. In fact, in the large LjEy region, the 
oscillations are averaged out and the survival probability there can be well approximated 
with 0.5 (for maximal mixing). In the region of small LjE^ both probabilities approach 
unity. In the region LjEy around 400 km/GeV, where the probability for the neutrino 
oscillation model has the first minimum, the two curves are most easily distinguishable, at 
least in principle. 

Decay Model 



There are two decay possibilities that can be considered: (a) z/2 decays to z/4 which is 
dominantly Ug with U2 and u^ mixtures of z/^ and u^-, as in Eq. ([T2|) , and (b) 1/2 decays into 
i?4 which is dominantly z/,- and z/2 and u^ are mixtures of z/^ and Ug. In both cases the decay 
interaction has to be of the form 

jO-int = 924 Ti^ iy2L-J + h.c. (15) 

where J is a Majoron field that is dominantly iso-singlet (this avoids any conflict with the 
invisible width of the Z). Viable models for both the above cases can be constructed pB| , pU| . 
However, case (b) needs additional iso-triplet light scalars which cause potential problems 
with Big Bang Nucleosynthesis (BBN), and there is some preliminary evidence from SuperK 
against z/^-z/^ mixing [^]. Hence we only consider case (a), i.e. z/2 — >■ z/4 + J with z/4 ^ Ug, 
as implicit in Eq. (p!2D. With this interaction, the z/2 rest-lifetime is given by 



IGvr 7712 

g"^ 5iv?{l + x) 



xu/i iib2 / \ 



where 8w? = m\ — m\ and x = 7714/1712 (0 < x < 1). From the value of a ^ = r2/?Ti2 = 
63 km/GeV found in the fit and for x = 0, we have 

/5m2~0.16eV^ (17) 

Combining this with the bound on g"^ from K ^ fi decays of g^ < 2.4 x 10^^ [^ we have 



6771' > 650 eV^ . (18) 

Even with a generous interpretation of the uncertainties in the fit, this 5m^ implies a min- 
imum mass difference in the range of about 25 eV. Then z/2 and z/s are nearly degenerate 
with masses > (9(25 eV) and z/4 is relatively light. We assume that a similar coupling of 
z/3 to z/4 and J is somewhat weaker leading to a significantly longer lifetime for z/3, and the 
instability of z/3 is irrelevant for the analysis of the atmospheric neutrino data. 

For the atmospheric neutrinos in SuperK, two kinds of tests have been proposed to 
distinguish between z/^-z/r oscillations and z/^-z/^ oscillations. One is based on the fact 



that matter effects are present for z/^-z/^ oscillations |^ but are nearly absent for z/^-z/. 



oscillations ^2[ leading to differences in the zenith angle distributions due to matter effects 



on upgoing neutrinos |]33|. The other is the fact that the neutral current rate will be affected 



in z/„-z/s oscillations but not for z/„-z/^ oscillations as can be measured in events with single 



33] . In these tests our decay scenario will behave as a hybrid in that there is no matter 



TfO's 

effect but there is some effect in neutral current rates. 

Long-Baseline Experiments 

The survival probability of z/^ as a function of L/E is given in Eq. (1). The conversion 
probability into u^ is given by 



This result differs from 
Furthermore, -P(z/^ 



P(^^ - l^r) 


= sin2^cos2 0(l-e-"^/2i?)2_ ^^g^ 


m 1 - P(z/^ 


-^ z/^) and hence is different from z/^-z/^ oscillations. 
> z/^) is not 1 but is given by 



Pil^^. ^ ^^,J + P(^^M ^ ^^r) = 1 - cos^ 9{1 - e-"^/^) (20) 

and determines the amount by which the predicted neutral-current rates are affected com- 
pared to the no oscillations (or the v^-Vt oscillations) case. In Fig. 3 we give the results for 
Pi^u^ -^ u^), Pi^u^ -^ Ut-) and Pi^u^ — > u^) + P{j-'^ -^ ^t) for the decay model and compare 
them to the Vp-'^^T oscillations, for both the K2K |^3[ and MINOS [^] (or the corresponding 



European project [^) long-baseline experiments, with the oscillation and decay parameters 
as determined in the fits above. 

The K2K experiment, already underway, has a low energy beam Ey ^ 1-2 GeV and a 
baseline L = 250 km. The MINOS experiment will have 3 different beams, with average 
energies E^ = 3, 6 and 12 GeV and a baseline L = 732 km. The approximate L/E^ ranges 
are thus 125-250 km/GeV for K2K and 50-250 km/GeV for MINOS. The comparisons in 
Figure 3 show that the energy dependence of u^ survival probability and the neutral current 
rate can both distinguish between the decay and the oscillation models. MINOS and the 
European project may also have r detection capabilities that would allow additional tests. 

Big Bang Nucleosynthesis 

The decay of 1^2 is sufficiently fast that all the neutrinos [ue, i^f^, Vri ^s) and the Majoron 
may be expected to equilibrate in the early universe before the primordial neutrinos decouple. 
When they achieve thermal equilibrium each Majorana neutrino contributes N^ = 1 and the 
Majoron contributes N^ = 4/7 [^], giving and effective number of light neutrinos N^ = 4| 
at the time of Big Bang Nucleosynthesis. From the observed primordial abundances of ^He 
and ^Li, upper limits on Ny are inferred, but these depend on which data are used p9|-^ 



Conservatively, the upper limit to N^ could extend up to 5.3 (or even to 6 if ^Li is depleted 
in halo stars p9[] ). 

Cosmic Neutrino Fluxes 

Since we expect both z/2 and 2/3 to decay, neutrino beams from distant sources (such as 
Supernovae, active galactic nuclei and gamma-ray bursters) should contain only Ue and z/g 
but no Up, Up, Vr and Vr- This is a very strong prediction of our decay scenario. We can 
compare the very different expectations for neutrino flavor mixes from very distant sources 
such as AGN's or GRB's. Let us suppose that at the source the flux ratios are typical 
of a beam dump, a reasonable assumption: N^^ : Ny : N^^ = 1:2:0. Then, for the 
conventional oscillation scenario, when all the 5m?^s satisfy 5m? L/4E » 1, it turns out 
curiously enough that for a wide variety of choices of neutrino mixing matrices, the final 
flavor mix is the same, namely: N^,^ : N^ : N^^ = 1:1:1. In the case of the decay B 
scenario, as mentioned here, we have N^^ : N^^ : N^^ = 1:0:0. The two are quite distinct. 
Techniques for determining these flavor mixes in future KM3 neutrino telescopes have been 
proposed |^ . 



Reactor and Accelerator Limits 

The Ue is essentially decoupled from the decay state z/2 so the null observations from the 
CHOOZ reactor are satisfied |^3|. The mixings of Up and Ur with z/^ and z/e are very small, so 



there is no conflict with stringent accelerator limits on flavor oscillations with large Sm [^ 

In summary, neutrino decay remains a viable alternative to neutrino oscillations as an 
explanation of the atmospheric neutrino anomaly. The model consists of two nearly de- 
generate mass eigenstates z/2, z/3 with mass separation > (9(25 eV) from another nearly 

8 



degenerate pair z/i, z/4. The z/^ and Vr flavors are approximately composed of z/2 and z/3, with 
a mixing angle 62^ — 57°. The state V2 is unstable, decaying to z/4 and a Majoron with a 
lifetime r2 ~ 10~^^ sec. The electron neutrino z/g and a sterile neutrino z/^ have negligible 
mixing with v^^v^- and are approximate mass eigenstates (z/e ~ viii's ~ ^4)5 with a small 
mixing angle 6'i4 and a 5m|]^ ^ 10^^ eV^ to explain the solar neutrino anomaly. The states z/3 
and z/4 are also unstable, but with z/3 lifetime somewhat longer and z/4 lifetime much longer 
than the z/2 lifetime. This decay scenario is difficult to distinguish from oscillations because 
of the smearing in both L and E^ in atmospheric neutrino events. However, long-baseline 
experiments, where L is fixed, should be able to establish whether the dependence of L/ E^ 
is exponential or sinusoidal. In our scenario only z/i is stable. Thus, neutrinos of supernovae 
or of extra galactic origin would be almost entirely z/g. The contribution of the electron neu- 
trinos and the Majorons to the cosmological energy density VL is negligible and not relevant 
for large scale structure formation. 

Another proposal for explaining the atmospheric neutrinos is based on decoherence of 
the z/'s in the fiux [^. The idea is that u' s are interacting and getting tagged before 



their arrival at the detector. The cause is unknown but could be a number of speculative 
possibilities such as a large neutrino background, new fiavor sensitive interactions in an extra 
dimension, violation of quantum mechanics etc. The z/^ survival probability goes as: 

p = [1 + cos2eexp{-t/T)] (21) 

The Super-Kamiokande data can be fit by choosing r ~ 10~^s and sin26 ~ 0.4. A detailed 
fit to all the data over the whole energy range has not been attempted yet. 



VII. CONCLUSION 

As I mentioned at the beginning, the main motivation for this exercise is to try to 
establish neutrino oscillations (due to mass-mixing) as the unique explanation of the observed 
anomalies. Even if neutrinos have masses and do mix, the observed neutrino anomalies may 
not be due to oscillations but due to other exotic new physics. These possibilities are testable 
and should be ruled out by experiments. I have tried to show that we are beginning to carry 
this program out. 
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FIG. 1. Comparison of decay model (solid histograms) and v^-Vr oscillation model (dashed 
histograms) with SuperK data. 
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FIG. 2. Survival probabiliity for i/^ versus logiQ{L/E) for the decay model (heavy solid curve) 
and ly^ oscillation model (thin curve). 
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FIG. 3. Long-baseline expectations for the K2K and MINOS long-baseline experiments from 
the decay model and the I'^-i^r oscillation model. The upper panel gives the neutral current 
predictions compared to no oscillations (or v^-Vr oscillations). 
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